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摘要 :在 全 球 变 暖 的 背景 下 ,由 土地 利用 变化 导致 的 土壤 碳 库 的 变化 已 经 受到 越 来 越 多 的 关注 。 首 先 采 用 物种 分 布 模型 预测 了 
天 山北 坡 雪 岭 云 杉林 的 潜在 分 布 ,其 次 估计 了 与 被 开明 为 农田 的 雪 岭 云 杉林 面积 (PSC) 以 及 由 林地 开 暴 为 农田 造成 的 有 机 碳 
损失 。PSC 分 别 由 雪 岭 云 杉 林 的 现实 分 布 潜在 分 布 和 农田 的 现实 分 布 确定 。 云 杉林 地 和 农田 的 土壤 有 机 碳 含量 由 野外 采样 
和 实验 室 分 析 获 得 。 研 究 发 现 ,PSC 面积 为 2.68x10"hm ,被 开垦 为 农田 的 雪 岭 云 杉 林 土 壤 有 机 碳 的 损失 为 171.7 v/hm? ;研究 
区 总 有 机 碳 的 损失 为 459.70Tg。 结 果 表明 ,研究 区 的 林地 恢复 和 重建 项 目 将 会 使 土壤 有 机 碳 储量 有 所 增加 , 且 土 壤 表 层 的 增加 
量 多 于 深层 。 

关键 词 :天 山 ; 雪 岭 云 杉林 ;潜在 分 布 ;土地 利用 变化 ;土壤 有 机 碳 损失 
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Estimation of soil organic carbon loss from Picea schrenkiana forest to farmland 
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Abstract: The effect of land-use changes on soil carbon stocks is an increasing concern in the context of global climate 
warming. In this study, we initially delineated the deforestation area of a Picea schrenkiana forest cultivated to cropland 
( PSC) at the north slope of the Tianshan Mountains. We then estimated the corresponding soil organic carbon (SOC) loss 
during the land-use change. The PSC was estimated by the potential, actual distributions of forest and cropland. The SOC 
content of the forest and cropland soils was obtained through field sampling and laboratory analysis. We found that the area 
of the PSC was 2.68 x 10^ hm^ , and the SOC loss ( per unit area) caused by the cultivation of forest to cropland was 171.7 
t/hm^, correspondingly. The total SOC loss in the study area was 459.70 Tg. This result implies that continuing the 
afforestation and reforestation programs currently implemented in the study area would increase the SOC. In addition, we 
found that the potential SOC sequestration is expected to increase to a greater extent in the upper soil layers than the lower 


layers. 
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在 过 去 的 几 十 年 里 ,大 气 中 温室 气体 的 浓度 与 全 球 变 暧 之 间 的 关系 一 直 是 科研 所 关注 的 焦点 问题 。 
Mauna Loa Observatory 和 Le Quere 的 观测 记录 表明 ,大 气 中 二 氧化 碳 ( C0, ,最 重要 的 温室 气体 ) 的 浓度 已 达到 
400ul/ L, 几乎 是 工业 革命 前 浓度 的 135% (278kLZL) 1。 联合 国政 府 间 气候 变化 专门 委员 会 (IPCC) 的 报 
告 表明 , 现 阶 段 大 气 中 CO, 浓度 和 相关 辐射 强迫 正在 以 前 所 未 有 的 速度 增加 中 ,上 且 大 气 中 C0, 浓 度 的 增加 主 
要 是 由 人 为 的 土地 利用 变化 和 化 石 燃料 的 使 用 造成 的 申 。 由 于 森林 土壤 中 储存 了 大 量 的 有 机 碳 ,因此 ,科学 
家 在 不 同 的 森林 生态 系统 中 ,开展 了 关于 森林 土地 利用 覆盖 变化 (LUCC) 和 相应 的 有 机 碳 损失 的 研究 |。 
然而 在 我 国 西 北 干旱 区 ,对 山地 森林 区 域 LUCC 造成 的 碳 损失 的 研究 较为 缺乏 。 

本 研究 估算 了 天 山北 坡 雪 岭 云 杉林 从 林地 到 耕地 的 土地 覆盖 变化 过 程 中 碳 的 损失 。 相 关 工 作 主 要 集中 
于 以 下 几 点 : (1) 对 遥感 图 像 进行 分 类 ,得 到 研究 区 的 雪 岭 云 杉林 和 农田 的 实际 分 布 ; (2) RH gc 
( Maximum entropy model，Maxent) 模 型 预测 雪上 岭 云 杉林 的 潜在 分 布 ;(3) 计 算 已 被 开 暴 为 农田 的 潜在 林地 的 
土地 面积 ( Potential Schrenk's spruce forest that has been cultivated to Cropland, PSC);(4) 采 集 土壤 样本 ,测定 
林地 和 耕地 中 土壤 有 机 碳 的 含量 和 密度 ;(5) 利 用 土壤 有 机 碳 的 含量 和 密度 数据 以 及 PSC 面积 估算 相关 土地 
利用 覆盖 变化 的 矶 损失 。 


1 材料 与 方法 


1.1 研究 区 域 

天 山 位 于 塔克拉玛干 沙漠 的 北部 和 西部 地 区 ,塔里木 岔 地 以 北 ( 图 1) ,是 喜马拉雅 造山 带 的 一 部 分 ,于 新 
生 代 时 期 由 印度 板块 和 亚 欧 板块 撞击 形成 。 天 山 的 植被 覆盖 类 型 (从 低 海 拔 到 高 海拔 ) 包括 草原 、 森 林 草 
原 、 亚 高 山 灌 从 、 高 山 草 馈 、 永 久 性 冰雪 。 天 山北 坡 ( 阴 坡 ) 草原 森林 带 的 范围 是 从 海拔 1400m 到 2800m ,和 森林 
主要 由 雪 岭 云 杉 ( Picea. schrenkiana) 林 组 成 。 


图 1 研究 区 示意 图 (天 山北 坡 , 蓝 色 的 阴影 区 域 ) 


Fig.l Study area (north slope of the Tianshan Mountains, blue-shaded regions) 


1.3 实际 林地 

研究 区 内 雪 岭 云 杉 林 和 农田 的 实际 分 布 是 通过 对 TM 3824/8 (2014 4E 6 H 1 H $8 H 31 日) 进行 分 类 
(分 类 和 回归 树 方 法 ,ENVI 4.8) 得 到 的 。 分 类 结果 的 精度 利用 214 个 随机 样 点 和 kappa 系数 进行 评价 ,通过 
计算 得 到 kappa 系数 为 0.8 ,表明 了 分 类 结果 的 合理 性 。 

L3 雪 岭 云 杉 林 的 潜在 分 布 模型 

本 研究 确定 了 155 个 雪 岭 云 杉 林 样 点 (139 个 位 于 中 国境 内 ,16 个 位 于 哈萨克 斯 坦 和 吉尔 吉 斯 斯 坦 ) ,用 
来 模拟 森林 物种 的 潜在 分 布 (图 1) 。 位 于 中 国 的 样 点 通过 2012—2014 年 的 野外 调查 确立 的 。 哈 萨克斯 坦 和 
吉尔 吉 斯 斯 坦 的 样 点 通过 谷歌 地 球 软件 确定 。 

雪 岭 云 杉林 的 潜在 分 布 由 Maxent 模型 进行 预测 ' ”1 。 该 模型 已 经 被 应 用 于 森林 的 恢复 规划 .气候 变化 
影响 预测 中等。 在 本 研究 中 ,选取 了 年 平均 气温 、 最 冷 月 份 的 最 低 气 温 、 年 平均 降水 量 、 最 湿润 季节 的 降水 
量 4 个 与 生物 气候 变量 以 及 由 ( Digital Elevation Model, DEM ) 计算 得 到 的 坡 向 ,坡度 ,太阳 辐射 量 和 地 形 湿度 
TEC P EON PRU E, RACES (3.3.3 K. 版 本 ) 在 运行 时 默认 值 的 设置 如 下 :随机 测试 百分比 ,25; 正 则 化 
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因子 ,1; 最 大 背景 样本 数 ,10000; 重 复 次 数 ,10。 实 际 运算 过 程 中 ,本 研究 使 用 了 75% 的 样 点 作为 训练 样本 ,其 
他 25% 为 验证 样本 。 获 得 雪 岭 云 杉林 潜在 分 布 的 概率 图 后 ,通过 最 大 化 敏感 性 和 特异 性 的 阔 值 选择 策略 将 
概率 值 转化 为 二 进 制 值 5m 。 

在 默认 的 情况 下 ,最 大 焙 模 型 采用 受 试 者 工作 特征 曲线 (Receiver Operating Character, ROC) 和 ROC 曲线 
下 面积 (Area Under ROC, AUC) 值 来 评估 模型 性 能 。 鉴 于 有 研究 对 AUC 提出 质疑 '” ,因此 ,在 本 研究 中 评价 
最 大 彤 模型 性 能 的 是 真实 技巧 统计 值 (True Skill Statistic, TSS) 。 根 据 Liu 等 的 研究 059 ,TSS 的 值 可 以 从 预测 
误差 矩阵 中 计算 得 到 。 
14 PSC 面积 估算 

被 开 晨 为 农田 的 潜在 林地 的 土地 面积 (PSC) 分 别 由 雪 岭 云 杉林 的 现实 分 布 (Actual distribution of the 
P.schrenkiana forests, Pa) 潜在 分 布 (Potential distribution of the P.schrenkiana forests , Pp ) 和 农田 的 现实 分 布 
( Actual distribution of the Cropland , Ca) 确定 ,可 表示 为 :Pae Pp, PpN Ca=PSC; H ArcGIS 计算 可 得 被 开 暴 为 
农田 的 潜在 林地 的 土地 面积 。 
1.5 URRIA 

野外 工作 中 ,在 每 一 个 野外 样 地 采样 点 (雪上 岭 云 杉林 和 农田 ) 平 行 的 采集 3 个 土壤 剖面 的 土壤 , 带 回 实 验 
室 分 析 获 得 该 样 点 土壤 有 机 碳 的 平均 值 。 在 每 个 剖面 以 10cm 为 间隔 采集 土壤 样本 (0 一 80cm ,或 遇 到 基 岩 为 
止 )。 实 验 室 分 析 过 程 中 ,采用 重 铬 酸 钾 容 量 法 测定 碳 的 密度 |。 获得 碳 含量 之 后 ,通过 如 下 公式 计算 得 到 
整个 土壤 剖面 的 有 机 碳 密度 5] ; 


NEL LU QUA (1) 


式 中 ,o 是 <2mm 土壤 颗粒 的 体积 密度 (g/cm ) ,8, 29 2 2mm 部 分 的 相对 体积 (%) , d, 为 i 土 层 的 厚度 (此 处 
d;=10cm) ,C; 是 i 土 层 的 兢 含 量 ( g/kg)。 


2 结果 


21 雪 岭 云 杉 林 实际 和 潜在 分 布 

雪 岭 云 杉林 的 实际 分 布 从 遥感 影像 的 分 类 结果 中 提取 。 结 果 表 明 , 天 山北 坡 雪 岭 云 杉林 的 面积 为 2.43x 
10^ hm’, Maxent 模型 模拟 结果 表明 雪 岭 云 杉 林 的 潜在 分 布 土地 面积 占 地 为 8.42x10* hm*(TSS 值 为 0.753, 表 
明 对 目标 物种 的 潜在 分 布 区 的 预测 是 合理 的 )。 图 2 显示 了 雪 岭 云 杉林 的 实际 分 布 和 潜在 分 布 。 通 过 比较 


图 2 雪 岭 云 杉林 和 实际 分 布 和 潜在 分 布 
Fig.2 Actual and pontential distributions the P.schrenkiana forests 
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雪 岭 云 杉林 的 实际 和 潜在 的 分 布 ,研究 发 现实 际 分 布 和 潜在 分 布 的 面积 差异 为 5.99x10" hm? ,实际 分 布 
(2.43x10* hm2) 仅 占 潜在 分 布 (8.42x105 hm?) 的 28.81%。 退 化 的 森林 可 能 已 经 被 开 星 为 农田 .草地 或 其 他 类 
型 的 土地 覆盖 或 土地 利用 类 型 。 本 研究 关注 的 是 雪 岭 云 杉 林 被 开垦 为 农田 潜在 分 布 。 经 过 计算 可 得 ,PSC 面 
FUA 2.68x10* hm^( KI 2) , 
2.2 雪上 岭 云 杉林 与 农田 的 土壤 有 机 碳 

室内 分 析 显 示 , 在 本 研究 的 采样 点 区 域 ,雪上 岭 云 杉林 表层 土壤 (0 一 10cm ) 的 土壤 有 机 碳 密 度 为 4.25 一 
13.79kg/m (42.5 一 137.9 ho ) ,平均 值 为 8.86kg/m^ (88.6 hm )( 表 1) 。 表 层 土壤 中 的 土壤 有 机 碳 密度 最 
高 ,之 后 土壤 有 机 碳 密度 沿 土壤 剖面 从 表层 至 底层 逐渐 降低 。 土 壤 放 面 深度 为 70 一 80cm 处 的 土壤 有 机 碳 密 
度 最 低 ,在 0.14—3.68kg/m" 之 间 (1.4 一 36.8 t/hm ) ,平均 为 1.88kg/m (18.8 V/hm^) 。 和 森林 土壤 有 机 碳 密度 为 
34.91kg/m/(349.1 t/hm?) 。 


表 1 雪 岭 云 杉林 土壤 剖面 的 土壤 有 机 碳 密度 
Table 1 SOC density of the soil profile of the P.schrenkiana forest 


m "em mom EX "m norm FHE 
T 土壤 深度 土壤 深度 à . 

i . SOC density/ Average/ SOC density/ Average/ 

: Soil depth/cm 2 A Soil depth/cm n A 
(kg/m") (kg/m) (kg/m) (kg/m") 

0—10 4.25—13.79 8.86 10—20 4.40—10.71 6.27 

20—30 2.05—10.73 4.87 30—40 1.69—8.24 3.47 

40—50 2.03—5.48 3.28 50—60 1.94—9.44 3.71 

60—70 0.19—4.47 2.57 70—80 0.14—3.68 1.88 


在 0 一 20cm 深度 的 农田 土壤 层 ,土壤 有 机 碳 密 度 6 
= (0 一 10cm JE X 2.40kg/m? , 10 一 20cm 土壤 层 为 
E 3.07kg/m , $ 1] Jy 2. 74kg/m? ) HE 20 一 50cm 土壤 层 
13 (20 一 30cm ,30—40cm ,40—50cm 土壤 层 分 别 为 1.66kg/ 
m’ ,2.40kg/m ,1.85kg/m ;平均 为 1.97kg/m ) 高 。 这 一 
结果 揭示 了 土壤 有 机 碳 密 度 从 表层 向 深层 的 递减 趋势 。 
M9 然而 ,减少 的 趋势 没有 延续 到 更 深 的 土壤 层 。 在 50— 
60cm 与 60—70cm 土 层 深度 ,土壤 有 机 碳 密度 各 自 增 加 
到 3.33 kg/m? 5j 3.02 kg/m? ,平均 密度 为 3.18 kg/m*( 图 


上 


土壤 有 机 碳 密度 
SOC density/(kg/n?) 
N w 


— 


3) 。 农 田 土壤 总 有 机 碳 密度 为 17.74 kg/m^ (177.4 t/ ”010 10—20 20—30 30—40 40—50 50—60 60—70 
2 土壤 深度 Soil depth/cm 

hm ) 。 

2.8. We 图 3 农田 土壤 齐 面 深度 与 土壤 有 机 碳 密度 之 间 的 关系 


由 前 述 结果 可 知 , 云 杉林 土壤 总 有 机 碳 密 度 为 Fig.3 General trend of the SOC density along the soil profiles of 
349.1 WVhm ,农田 土壤 总 有 机 碳 密 度 为 177.4 t/hm”。 ‘Pe cropland 
因此 ,由 林地 到 农田 的 土地 利用 覆盖 变化 ,可 导致 土壤 
总 有 机 碳 含量 的 变化 量 为 171.7 thm , 约 为 林地 土壤 总 有 机 碳 的 近 50%( 49.18%)。 当 在 整个 研究 区 内 考虑 
土壤 的 碳 损失 时 ( 即 考虑 PSC 土地 面积 2.68x 10* hm?) ,计算 可 得 从 森林 土壤 到 农田 对 应 的 土地 覆盖 变化 的 
碳 损失 约 为 459.70Tg。 然 而 ,需要 注意 的 是 ,林地 转化 为 农田 之 后 ,与 自然 土壤 相 比 ,农田 土壤 碳 密度 的 变化 
是 自然 因素 和 农业 管理 的 综合 作用 。 考 虑 到 本 研究 中 采集 的 土壤 样本 来 自 于 农田 经 营 管理 措施 之 下 ,相关 研 
究 表明 通过 秸秆 还 田 与 有 机 肥 施 用 、 化 肥 投入 增加 与 合理 的 养分 配 比 以 及 少 ( 免 ) 耕 技术 的 推广 ,会 导致 农田 
土壤 有 机 碳 的 增加 。 因 此 ,本 文 计算 的 林地 与 农田 总 有 机 碳 密 度 的 差异 可 能 会 减 小 ,也 即 最 终 的 碳 损失 量 
可 能 会 有 所 减少 。 
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由 于 各 种 环境 和 人 为 因素 的 影响 ,难以 获得 林地 被 
耕作 后 的 土壤 有 机 碳 损失 量 。Niu 和 Duiker 通过 在 林 3 A. 
地 边缘 农业 用 地 上 划 定 潜在 林地 来 评定 实践 造林 的 碳 j 
损失 。 在 估算 过 程 中 ,作者 利用 土地 歼 盖 数 据 集 和 m | y 
土壤 地 理 数据 库 来 确定 潜在 的 林业 用 地 。 这 不 失 为 一 E31 
种 估计 碳 损 失 量 的 方法 ,然而 , 划 定 的 林业 用 地 在 生物 ET | 
气候 空间 上 是 否 适 宜 森林 生长 仍 不 清楚 。Schulp 5 等 R2 ,| 
将 整个 欧洲 作为 研究 区 开展 了 未 来 土地 利用 变化 对 土 ““ ,| 
壤 和 植被 碳 储量 变化 影响 的 研究 。 在 该 研究 中 ,作者 将 | 7 
土地 利用 变化 模型 (Dyna-CLUE ) 和 碳 模 型 (大 尺度 l l l 
bookkeeping 模型 ) 相 结合 来 预测 未 来 的 土地 利用 和 相 RE Co E NE 
c» ”应 的 碳 损失 量 。 该 研究 对 比 了 不 同 尺度 上 两 种 模型 进 
行 结合 所 关联 的 复杂 性 。 鉴 于 土地 利用 变化 模型 和 土 图 4 森林 和 农田 土壤 有 机 碳 的 差异 


地 利 用 情景 预测 的 不 确 定性 ， 造 成 评估 结 果 的 变化 较 Fig. 4 Difference in SOC density between the forest and 
大 。 在 本 研究 中 ,通过 实地 调查 实验 室 分 析 和 模型 模 mat soi 
拟 估 算 了 天 山北 坡 农田 的 碳 损失 量 。 与 基于 过 程 估算 
相 比 ,本 研究 所 要 求 的 数据 类 型 较 少 ,因此 适用 于 在 区 域 范 围 内 做 相应 的 估算 。 然 而 需要 指出 的 是 ,由 于 
SDM 模型 具有 不 确定 性 ,所 以 ,提高 潜在 分 布 模拟 精度 也 是 必要 之 举 。 

森林 和 农田 土壤 有 机 碳 的 差异 为 可 能 存在 的 碳 汇 提供 了 科学 依据 。 在 土地 利用 变化 及 其 土壤 有 机 碳 储 
量变 化 的 研究 中 ,Don ”等 认为 最 高 的 碳 损 失 是 由 原始 森林 转变 为 农田 ( 2596) 和 多 年 生 农作物 ( 3006 ) 3 
成 的 。 本 研究 结果 发 现 如 果农 田 能 够 恢复 成 天 然 林 ,那么 可 能 带 来 矶 储量 的 增加 。 由 于 土壤 有 机 碳 的 增加 可 
能 发 生 在 每 一 土 层 中 ,因此 研究 土壤 有 机 碳 在 土壤 各 层 中 的 增加 差异 ,有 助 于 深入 理解 土地 利用 变化 对 应 的 
土壤 有 机 碳 的 变化 过 程 。 本 研究 除了 进一步 证 实 林地 土壤 有 机 碳 随 着 土壤 深度 的 增加 而 降低 的 趋势 ,还 
发 现 森 林 与 农田 上 层 土壤 有 机 碳 含量 的 差异 大 于 下 层 土壤 (图 4) ,因此 上 层 土壤 可 能 会 更 多 的 贡献 于 退耕 还 
O ” 林 及 其 他 森林 恢复 工程 所 产生 的 碳 汇 增加 。 
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